ON THE DISTRIBUTIONAL HESSIAN OF THE DISTANCE FUNCTION 
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ABSTRACT. We describe the precise structure of the distributional Hessian of the distance function 
from a point of a Riemannian manifold. In doing this we also discuss some geometrical properties 
of the cutlocus of a point and we compare some different weak notions of Hessian and Laplacian. 



1. Introduction 

Let (M, g) be an n-dimensional, smooth, complete Riemannian manifold, for any point p € M 
we define d p : M — > M to be the distance function from p. 

Such distance functions and their relatives, the Busemann functions, enters in several argu- 
ments of differential geometry. It is easy to see that, apart from the obvious singularity at the 
point p, with some few exceptions such distance function is not smooth in M \ {p} (for instance, 
when the manifold M is compact), it is anyway 1-Lipschitz and differentiable with a unit gradi- 
ent almost everywhere (by Rademacher's theorem). 

In this note we are concerned with the precise description of the distributional Hessian of d p , 
having in mind the following Laplacian and Hessian comparison theorems (see [25], for instance). 

Theorem 1.1. If (M, g) satisfies Ric > (n — 1)K then, considering polar coordinates around the points 
p G M and P in the simply connected, n-dimensional space S K of constant curvature K G R, we have 

Ad p (r) < A K d^(r). 

If the sectional curvature of '(A/, g) is greater or equal to K, then 

Hess d p (r) < Hcss^df (r) . 

Here A K dp(r) and Hcss K dp (r) denote respectively the Laplacian and the Hessian of the distance func- 
tion d'^(-) = d K (P, ■) in S K , at distance r from P. 

It is often stated by several authors that these inequalities actually hold on the whole mani- 
fold (M, g), in some weak sense, that is, in sense of distributions, of viscosity, of barriers. Such 
conclusion can simplify and sometimes is actually necessary in global arguments involving this 
comparison theorem, more in general, one often would like to use (weak or strong) maximum 
principle for the Laplacian in situations where the functions involved are not smooth, for instance, 
in the proof of the "splitting" theorem (first proved by Cheeger and Gromoll [9]) by Eschenburg 
and Heintze [13], but also of Topogonov theorem and the "soul" theorem (Cheeger, Gromoll and 
Meyer [10, 17]). 

To be precise, we give the respective definitions of these notions. 

Definition 1.2. Let A be a smooth, symmetric 2-form on a Riemannian manifold (Af, g). 

• We say that the function / : M — >• R satisfies Hess / < A in distributional sense if for 
every smooth vector field V with compact support there holds f M f V^T^V^) c?Vol < 
f M A ij V i Vi<No\. 

• For a continuous function / : AI — > R, we say that Hess / < A at the point p G M in barrier 
sense if for every e > there exists a neighborhood U £ of the point p and a C 2 -function 
h £ : U £ — > R such that h e (p) = f(p), h e > f in U £ and Hess h £ (p) < A(p)+eg(p) as 2-forms 
(such a function h £ is called an upper barrier). 
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• For a continuous function / : M — >• R, we say that Hess / < A at the point p e M in 
viscosity sense if for every C 2 -function ft, from a neighborhood t/ of the point p such that 
ft(p) = f(p) and ft < / in U, we have Hess h(p) < A(p). 

The weak notions of the inequality A/ < a, for some smooth function a : M — > R, are defined 
analogously. 

• We say that the function / : M — > R satisfies A/ < a in distributional sense if for 
every smooth, nonnegative function y> : ill — > R with compact support there holds 
J M fA<pdVol<f M o^dVoL 

• For a continuous function / : M — > R, we say that A/ < a at the point p € M in barrier 
sense if for every e > there exists a neighborhood U e of the point p and a C 2 -function 
ft £ : U e -> R such that h £ (p) = f(p), h £ > f in J7 e and Aft £ (p) < a(p) + e. 

• For a continuous function / : M — > R, we say that A/ < a at the point p G M in 
viscosity sense if for every C 2 -function ft from a neighborhood £/ of the point p such that 
ft(p) = f(p) and ft < / in U, we have Aft(p) < a(p). 

The notion of inequality "in barrier sense" was defined by Calabi [8] (for the Laplacian) back 
in 1958 (he used the terminology "weak sense" rather than "barrier sense") who also proved 
the relative global "weak" Laplacian comparison theorem (see also the book of Petersen [25, Sec- 
tion 9.3]). 

The notion of viscosity solution (which is connected to the definition of inequality "in viscosity 
sense", see Appendix A) was introduced by Crandall and Lions [12, Definition 3.2] for partial 
differential equations, the above definition for the Hessian is a generalization to a very special 
system of PDEs. 

The distributional notion is useful when integrations (by parts) are involved, the other two con- 
cepts when the arguments are based on maximum principle. 

It is easy to see, by looking at the definitions, that "barrier sense" implies "viscosity sense", more- 
over, by the work [19], if / : M — > R satisfies A/ < a in viscosity sense it also satisfies Af < a 
as distributions and viceversa. In the Appendix A we will discuss in detail the relations between 
these definitions. 

In the next section we will describe the distributional structure of the Hessian (and hence of 
the Laplacian) of d p which will imply the mentioned validity of the above inequalities on the 
whole manifold. 

It is a standard fact that the function d p is smooth in the set M \ ({p} U Cut p ), where Cut p is the 
cutlocus of the point p, which we are now going to define along with stating its basic properties 
(we keep the books [16] and [26] as general references). It is anyway well known that Cut p is 
a closed set of zero (canonical) measure. Hence, in the open set M \ ({p} U Cut p ) the Hessian 
and Laplacian of d p are the usual ones (even seen as distributions or using other weak defini- 
tions) and all the analysis is concerned to what happens on Cut p (the situation at the point p is 
straightforward as d p is easily seen to behaves as the function ||x|| at the origin of R"). 

We let U p = {v € T P M | g p (v, v) = 1} to be the set of unit tangent vectors to M at p. Given 
v G Up we consider the geodesic 7„(£) = cxp p (tv) and we let a v e R + (possibly equal to +oo) 
to be the maximal time such that 7„([0, cr v ]) is minimal between any pair of its points. It is so 
defined a map a : U p —t R + U {+00} and the point j v (a v ) (when a v < +00) is called the cutpoint 
of the geodesic ~y v . 

Definition 1.3. The set of all cutpoints ~f v (a v ) for v e U p with a v < +00 is called the cutlocus of 
the point p £ M. 

The reasons why a geodesic ceases to be minimal are explained in the following proposition. 

Proposition 1.4. If for a geodesic j v (i) from the point p g M we have a v < +00, at least one of the 
following two (mutually non exclusive) conditions is satisfied: 

(1) at the cutpoint q = j v (a v v) there arrives at least another minimal geodesic from p, 

(2) the differential dexp p is not invertible at the point a v v <G T P M. 

Conversely, if at least one of these conditions is satisfied the geodesic j v (t) cannot be minimal on an interval 
larger that [0, a v ]. 
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It is well known that the subset of points q G Cut p where more than a minimal geodesic from 
p arrive coincides with Sing, which is the singular set of the distance function d p in M \ {p}. We 
also define the set Conj of the points q = j v (<t v ) G Cut p with e?exp p not invertible at a v v G T p M, 
we call Conj the locus of optimal conjugate points (see [16, 26]). 

Acknowledgments. We thank Giovanni Alberti, Luigi Ambrosio, Davide Lombardo and Federico Poloni 
for several valuable suggestions. The first author is partially supported by the Italian project FIRB-IDEAS 
"Analysis and Beyond". 

2. The Structure of the Distributional Hessian of the Distance Function 

The following properties of the function d p and of the cutlocus of p G M are proved in the 
paper [22], Section 3 (see also the wonderful work [20] for other fine properties, notably the local 
Lipschitzianity of the function a : U p — ► R + U {+oo}, in Theorem 1.1). 

Given an open set O c R", we say that a continuous function u : il — >• R is locally semiconcave 
if, for any open convex set K C SI with compact closure in fi, the function u\k is the sum of a C 2 
function with a concave function. 

A continuous function u : M —> R is called locally semiconcave if, for any local chart %jj : R" — >• U C 
M, the function u o ijj is locally semiconcave in R" according to the above definition. 

Proposition 2.1 (Proposition 3.4 in [22]). The function d p is locally semiconcave in AT \ {p}. 

This fact, which follows by recognizing d p as a viscosity solution of the eikonal equation | Vit| = 1 
(see [22]), has some relevant consequences, we need some definitions for the precise statements. 

Given a continuous function u : Q — > R and a point q G M, the superdifferential of u at q is the 
subset of T*M defined by 

d + u(q) = jcftp(g) | <p G C 1 (M),cp(q) - u(q) = min^ — u | • 

For any locally Lipschitz function u, the set d + u(q) is a compact convex set, almost everywhere 
coinciding with the differential of the function u, by Rademacher's theorem. 

Proposition 2.2 (Proposition 2.1 in [1]). Let the function u : M — > R be semiconcave, then the superdif- 
ferential d + u is not empty at each point, moreover, d + v is upper semicontinuous, namely 

qk -> q, v k ->• v, v k G d + u(q k ) =^> v G d + u(q). 
In particular, if the differential du exists at every point of M, then u G C 1 (M). 

Proposition 2.3 (Remark 3.6 in [1]). The set Ext{d + d p {q) of extremal points of the (convex) superdif- 
ferential set of d p at q is in one-to-one correspondence with the family Q{q) of minimal geodesies from p 
to q. Precisely Q (q) is described by 

Q(q) = { exp,(-trf), * G [0, 1] | V« G Ext(d+d p (q) } . 

We now deal with structure of the cutlocus of p G M. Let denote the (n— l)-dimensional 

Hausdorff measure on (M, g) (see [14, 27]). 

Definition 2.4. We say that a subset S C M is C r -rectifiable, with r > 1, if it can be covered by 
a countable family of embedded C -submanifolds of dimension (n — 1), with the exception of a 
set of H rl_1 -zero measure (see [14, 27] for a complete discussion of the notion of rectifiability). 

Proposition 2.5 (Theorem 4.10 in [22]). The cutlocus ofp G M is C°° -rectifiable. Hence, its Hausdorff 
dimension is at most n — 1. Moreover, for any compact subset K of M the measure - H n_1 (Cut p n K) is 
finite (Corollary 1.3 in [20]). 

To explain the following consequence of such rectifiability, we need to introduce briefly the 
theory of functions with bounded variation, see [5, 7, 14, 27] for details. We say that a function 
«:R"-> R m is a function with locally bounded variation, that is, u G BV\ ac , if its distributional de- 
rivative Du is a Radon measure. Such notion can be easily extended to maps between manifolds 
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using smooth local charts. 

A standard result says that the derivative of a locally semiconcave function stays in BV\ oc , in view 
of Proposition 2.1 this implies that the vector field Vd p belongs to BVi oc in the open set M \ {p}. 

Then, we define the subspace of BV\ oc of functions (or vector fields, as before) with locally 
special bounded variation, called SBV\ oc (see [2, 3, 4, 5, 7]). 

The Radon measure representing the distributional derivative Du of a function u : R" — > W n with 
locally bounded variation can be always uniquely separated in three mutually singular measures 

Du=- Du + Jit + Cu 

where the first term is the part absolutely continuous with respect to the Lebesgue measure C n , 
Ju is a measure concentrated on an (n — l)-rectifiable set and Cu (called the Cantor part) is a 
measure which does not charge the subsets of Hausdorff dimension (n — 1). 
The space SBVi oc is defined as the class of functions u G BV\ oc such that Cu = 0, that is, the 
Cantor part of the distributional derivative of u is zero. Again, by means of the local charts, this 
notion is easily generalized to Riemannian manifolds. 

Proposition 2.6 (Corollary 4.13 in [22]). The (W 1-1 -almost everywhere defined) measurable unit vector 
field Vd p belongs to the space SBV\ oc {M \ {p}) of vector fields with locally special bounded variation. 

The immediate consequence of this proposition is that the 2-f orm distribution Hess d p is actu- 
ally a Radon measure with an absolutely continuous part, with respect to the canonical volume 
measure Vol of (M, g), concentrated in M \ ({p} U Cut p ) where d p is a smooth function, hence 
in this set Hess d p coincides with the standard Hessian 2-form Hess d p times the volume measure 
Vol. When the dimension of M is at least two, the singular part of the measure Hess d p does not 
"see" the singular point p, hence, it is concentrated on Cut p , absolutely continuous with respect 
to the Hausdorff measure % n ~ l , restricted to Cut p . 

By the properties of rectifiable sets, at H n_1 -almost every point q G Cut p there exists an {n— 1)- 
dimensional approximate tangent space apT 9 Cut p c T q M (in the sense of geometric measure theory, 
see [14, 27] for details). To give an example, we say that anhyperplane T C K" is the approximate 
tangent space to an (n — 1) -dimensional rectifiable set K G M. n at the point xo, if W 1 ^ 1 L T is the 
limit in the sense of Radon measures, as p — > +oo, of the blow-up measures TL n ~ l L p(K — xq) 
around the point xq. With some technicalities, this notion can be extended also to Riemannian 
manifolds. 

Moreover, see [5], at % n_1 -almost every point q G Cut p , the field Ve? p has two distinct approximate 
(in the sense of Lebesgue differentiation theorem) limits "on the two sides" of apTgCut p c T q M, 
given by Vd+ and Vd p . 

We want to see now that at 'H n_1 -almost every point of Cut p there arrive exactly two distinct 
geodesies and no more. We underline that a stronger form of this theorem was already obtained 
in [6] and [15], concluding that the set Cut p \ U (where U is like in the following statement) 
actually has Hausdorff dimension not greater that n — 2. 

Theorem 2.7. There is an open set U c M such that •H"" 1 (Cut p \ U) = and 

• the subset Cut p n U does not contain conjugate points, hence the set of optimal conjugate points 
has W. n ~ 1 -zero measure; 

• at every point o/Cut p n U there arrive exactly two minimal geodesies from p G M; 

• locally around every point of Cut p n U the set Cut p is a smooth (n — l)-dimensional hypersurface, 
hence, apT g Cut p is actually the classical tangent space to a hypersurface. 

Proof. First we show that the set of optimal conjugate points Cory is a closed subset of W l ~ 1 - 
zero measure, then we will see that the points of Sing \ Conj where there arrive more than two 
geodesies is also a closed subset of "H n_1 -zero measure. The third point then follows by the 
analysis in the proof of Proposition 4.7 in [22]. 

Recalling that U p = {v G T p M \ g p (v, v) = 1} is the set of unit tangent vectors to M at p, we 
define the function c : U p — >• M+ U {+00} such that the point Y„(c-y) is the first conjugate point 
(if it exists) along the geodesic j v , that is, the differential dcxp p is not invertible at the point 
c v v G TpM. By Lemma 4.11 and the proof of Proposition 4.9 in [22], in the open subset V C U p 
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where the rank of the differential of the map F : U p — > M, defined as F(v) = exp (c v v) is n — 1, 
the map c : {/ p — > K + U {+00} is smooth hence .F(V) is locally a smooth hypersurface. As, by 
Sard's theorem, the image of U p \ V is a closed set of % n_1 — zero measure, we only have to deal 
with the images F(v) of the unit vectors v £ V with c„ = <7„ (see at the end of the introduction), 
that is, with F(V) n Cut p , which is a closed set. 

We then consider the set D c (F(V) n Cut p ) of the points q where apT g Cut p exists and the density 
of the rectifiable set F(V) n Cut p in the cutlocus of the point p, with respect to the Hausdorff 
measure 'H n ~ 1 , is one (see [14, 27]). It is well known that D and F(V) n Cut p only differ by a set 
of %" _1 -zero measure. If F(v) = q £ D then, c v = a v and, by the above density property, the 
hypersurface F(V) is "tangent" to Cut p at the point q, that is, T q F(V) = apT 9 Cut p . 
We claim now that the minimal geodesic j v is tangent to the hypersurface F(V), hence to the 
cutlocus, at the point q. Indeed, as dexp p is not invertible at c v v £ T P M, by Gauss lemma there 
exists a vector w £ T V U P such that dcxp p [c v v](w) = 0, hence 

dF v (w) = (dc[w](iy))7„(c v ) + d exp p [c v v](c v w) = (dc[i;](w))7„(c„) , 

thus, 7t)(c„) belongs to the tangent space dF(T v U p ) to the hypersurface F(V) at the point q, which 
coincides with apT^Cutp, as we claimed. 

By the properties of SBV functions described before, at % n -almost every point q £ D, the 
blow-up of the function d p is a "roof", that is, there arrive exactly two minimal geodesies both 
intersecting transversally the cutlocus at q (the vectors Vc£+ and Vd~ do not belong to &pT q M), 
hence the above minimal geodesic 7„ cannot coincide with any of these two. 
We then conclude that W l ^ 1 {D) = and the same for the set Conj. 

Suppose now that q £ Cut p \ Conj C Sing, by the analysis in the proof of Proposition 4.7 in [22] 
(and Lemma 4.8), at the point q there arrive a finite number m > 2 of distinct minimal geodesies 
and when m > 2 the cutlocus of p is given by the union of at least m smooth hypersurfaces with 
Lipschitz boundary passing at the point q, in particular the above blow-up at q cannot be a single 
hyperplane apT 9 Cut p . By the above discussion, such points with m > 2 are then of 'H"~ 1 -zero 
measure, moreover, by Propositions 2.2 and 2.3 the set of points in Cut p \ Conj with only two 
minimal geodesies is open and we are done. □ 

Remark 2.8. In the special two-dimensional and analytic case, it can be said something more, that 
is, the number of optimal conjugate points is locally finite and the cutlocus is a locally finite graph 
with smooth edges, see the classical papers by Myers [23, 24]. We conjecture that, in general, the 
set of optimal conjugate points is an (n — 2)-dimensional rectifiable set. 

By the third point if this theorem, in the open set U the two side limits Ve?+ and Vd~ of the 
gradient field Vd p are actually smooth and classical limits, moreover it is locally defined a smooth 
unit normal vector v q e T q M orthogonal to T«jCut p , with the convention that g q {v q ,v) is positive 
for every vector v £ T q M belonging to the half space corresponding to the side associated to Vd p . 
Hence, since "H n-1 (Cut p \ U) = 0, we have a precise description of the singular "jump" part as 
follows, 

JVrf p = -( (W+ - Vd p ) <g> i/) H 71 - 1 L Cut p 

and, noticing that the "jump" of the gradient of d p in U must be orthogonal to the tangent space 
T q Cut p , thus parallel to the unit normal vector v q £ T q M, hence we conclude 

JVrf p = -(i/«i/)|Vd+- Wd p \ g U n - 1 LCut p . 

Notice that the singular part of the distributional Hessian of d p is a rank one symmetric 2-form. 

Remark 2.9. This description of the "jump" part of the singular measure is actually a direct conse- 
quence of the structure theorem of BV functions (see [5]), even without knowing, by Theorem 2.7, 
that the cutlocus is 'H n_1 -almost everywhere smooth. 

Theorem 2.10. Ifn>2, the distributional Hessian of the distance from a point p £ M is given by the 
Radon measure 

Hess d p = Hess d p Vol -(v®v) | Vd+ - Vd p \ g U n ~ l L Cut p , 



6 



CARLO MANTEGAZZA, GIOVANNI MASCELLANL AND GENNADY URALTSEV 



where Hess d p is the standard Hessian 2-form ofd p , where it exists CH n 1 -almost everywhere on M), and 
Vdp, Vd~, v are defined above. 

Corollary 2.11. Ifn>2, the distributional Laplacian ofd p is the Radon measure 

Ad p = Ad p Vol - |Vrf+ - Vd" \ g H n ~ 1 L Cut p , 
where Ad p is the standard Laplacian of d p , where it exists. 
Corollary 2.12. There hold 

Ad p < Adp Vol 

and 

Hess dp < Hess d p Vol , 

as 2-forms. 

As a consequence, the Hessian and Laplacian inequalities in Theorem 1.1 hold in the sense of distributions. 
Moreover, we have 

Ad p > Adp Vol - 2U n ~ 1 L Cut p . 

and 

Hess dp > Hess d p Vol - 2{v ® v) H n ~ 1 L Cut p > Hess d p Vol - 2g "H"" 1 L Cut p , 
as 2-forms 

Remark 2.13. By their definition, it is easy to see that the same inequalities hold also for the Buse- 
mann functions, see for instance [25, Subsection 9.3.4] (in Section 9.3 of the same book it is shown 
that the above Laplacian comparison holds on the whole M in barrier sense while the analogous 
result for the Hessian can be found in Section 11.2). We underline here that Propositions 2.1, 2.2 
and 2.3 about the semiconcavity and the structure of the superdifferential of the distance function 
dp can also be used to show that the above inequalities hold in barrier /viscosity sense. 

Remark 2.14. Several of the conclusions of this paper holds also for the distance function from a 
closed subset of M with boundary of class C 3 at least, see [22] for details. 

Appendix A. Weak Definitions of Sub/Supersolutions of PDEs 

Let (M, g) be a smooth, complete, Puemannian manifold and let A be a smooth 2-form. 

If / : M — >■ R satisfies Hess / < A at the point p e M in barrier sense, for every e > there 
exists a neighborhood U e of the point p and a C 2 -function h E : U e — >• R such that h e {p) = f(p), 
h e > / m U e and Hess h £ (p) < A(p) + eg(p), hence, every C 2 -function h from a neighborhood U 
of the point p such that h(p) — f(p) and h < f in t/ satisfies h(p) = h e (p) and h < h £ in U fl U e . It 
is then easy to see that it must be Hess h(p) < Hess h e (p) < A(p) + eg(p), for every e > 0, hence 
Hess h(p) < A{p). This shows that Hess / < A at the point p e M also in viscosity sense. 
The converse is not true, indeed, it is straightforward to check that the function / : K — > R given 
by /(.t) = x 2 sin (1/x) when x ^ and /(0) = satisfies /"(0) < in viscosity sense but not in 
barrier sense. 

The same argument clearly also applies to the two definitions of Af < a, for a smooth function 

a : M ->• M. 

We see now that instead the definitions of viscosity and distributional sense coincide. 

Proposition A.l. If f : M — > R satisfies Hess / < A in viscosity sense, it also satisfies Hess / < A in 
distributional sense and viceversa. The same holds for Af < a. 

In order to show the proposition, we recall the definitions of viscosity (sub/ super) solution to 
a second order PDE. Take a continuous map F : ft x R x R" x S" —> R, where 51 is an open subset 
of R" and S n denotes the space of real nxn symmetric matrices; also suppose that F satisfies the 
monotonicity condition 

X>Y F(x,r,p,X)<F(x,r,p,Y), 

for every (x, r,p) <G O x R x R", where X > Y means that the difference matrix X — Y is 
nonnegative definite. We consider then the second order PDE given by F(x, f, V/, V 2 /) = 0. 
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A continuous function / : f2 — > R is said a viscosity subsolution of the above PDE if for every 
point xe Oand^G C 2 (fl) such that f(x) — <p(x) = sup n (/ — there holds F(x, <ys, V<ys, V 2 (p) < 
(see [11, 19]). Analogously, / e C°(Q) is a viscosity supersolution if for every point x <E and </? g 
C 2 {&) such that /(x)-y>(x) = inf n (/ -p), there holds F(x,(p,V(p,V 2 ip) > 0. If / e C°(fi) is both 
a viscosity subsolution and supersolution, it is then a viscosity solution of F(x, /, V/, V 2 /) = in 
ft. 

It is easy to see that the functions / e C°(fi) such that A/ < a in viscosity sense at any point of 
O, as in Definition 1.2, coincide with the viscosity supersolutions of the equation — A/ + a = at 
the same point (here the function F is given by F(x, r, p, X) = — trace X + a(x)). 

In the case of a Riemannian manifold (M, g), one works in local charts and the operators we 
are interested in become 

He</(x)=!^-r*,(x)|4 
A^(x)=.^(x)Hess£//(x), 
where T*- are the Christoffel symbols. 

Analogously to the case of R™, taking F(x, r,p,X) = -g ZJ (x)Xij + g lJ (x)T^ (x)pk + a(x) (which 
is a smooth function independent of the variable r), we see that, according to Definition 1.2, / 
satisfies A M / < a in viscosity sense at any point of M if and only if it is a viscosity supersolution 
of the equation F(x, /, V/, V 2 /) = at the same point. 

Getting back to R", given a linear, degenerate elliptic operator L with smooth coefficients, that 
is, defined by 

Lf(x) = -a i3 {x)Vlf{x) + b k (x)V k f(x) + c(x)f(x) , 

and a smooth function a : Q — > R, we say that / S C°(Q) is a distributional supersolution of the 
equation Lf + a = when 

(jL*ip + aip ^jdx > 

for every nonnegative, smooth function ip e C7^°(0). Here L* is the formal adjoint operator of L: 

L*y(x) = -V^a'VHx) - V k (b k cp)(x) + c(xMx). 

Under the hypothesis that the matrix of coefficients (a^ ) (which is nonnegative definite) has 
a "square root" matrix belonging to C 1 (f2, S n ), Ishii showed in paper [19] the equivalence of the 
class of continuous viscosity subsolutions and the class of continuous distributional subsolutions 
of the equation Lf + a = 0. More precisely, he proved the following two theorems (see also [21]). 

Theorem A.2 (Theorem 1 in [19]). If f e C°(Q) is a viscosity subsolution of the equation Lf + a = 0, 
then then it is a distribution subsolution of the same equation. 

Theorem A.3 (Theorem 2 in [19]). Assume that the "square root" of the matrix of coefficients (a,j) 
belongs to C 1 (il). If f e C°(f2) is a distributional subsolution of the equation Lf + a = 0, then then it 
is a viscosity subsolution of the same equation. 

As the PDE is linear, a function / e C°(£l) is a viscosity (distributional) supersolution of the 
equation Lf + a = if and only if the function — / is a viscosity (distributional) subsolution of 
L(—f) — a — 0, in the above theorems every occurrence of the term "subsolution" can replaced 
with "supersolution" (and actually also with "solution"). 

For simplicity, we will work in a single coordinate chart of M mapping onto SI C R™, while 
the general situation can be dealt with by means of standard partition of unity arguments. 
Consider / e C°(M) which is a viscosity supersolution of — A / + a = 0. It is a straightforward 
computation to check that this happens if and only if / is a viscosity supersolution of —yfg~A AI f + 
a^fg = 0, where ^fg = det g^- is the density of Riemannian volume of (M, g), and viceversa. 
Moreover, notice that setting L = —^/gA M we have that L* = L, that is, L is a self-adjoint 
operator; it also satisfies the hypotheses of Ishii's theorems, being the matrices and g lJ smooth 
and positive definite in O (see [18, Chapter 6], in particular Example 6.2.14, for instance). 
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Then, in local coordinates, Ishii's theorems guarantee that / is a distributional supersolution 
of the same equation, that is, / satisfies, for each cp g (SI), 



hence, 



/ fL*<pdx> — / cty/gipdx, 
Jn Jn 

g 

[ -fA M LpdVo\= I -f^gA M Lpdx> - I a^/g<pdx = - [ aipdVol. 
Jm Jn Jn Jm 



This shows that then / satisfies A M f < a in distributional sense, as in Definition 1.2. 

Following these steps in reverse order, one gets the converse. Hence, the notions of A M < a 
in viscosity and distributional sense coincide. 

Now we turn our attention to the Hessian inequality; it is not covered directly by Ishii's theo- 
rems, which are peculiar to PDEs and do not deal with systems (like the general theory of viscosity 
solutions). For simplicity, we discuss the case of an open set SI c M. n (with its canonical flat met- 
ric), since all the arguments can be extended to any Riemannian manifold (M, g) by localization 
and introduction of the first-order correction given by Christoffel symbols, as above. 

The idea is to transform the matrix inequality Hess f < A into a family of scalar inequalities; 
indeed, if everything is smooth, such inequality is satisfied if and only if for every compactly- 
supported, smooth vector field W we have W' l W : > HesSy / < AijW' L WK The only price to pay is 
that we lose the constant coefficients of the Hessian, hence making the linear operator L w , acting 
on / e C 2 (fl) as L w f — —W % Wi Hess^ /, only degenerate elliptic. Notice that Ishii's condition in 
Theorem A. 3 is satisfied for every smooth vector field W such that \\W\\ e C}. (SI), but not by any 
arbitrary smooth vector field. This has the collateral effect of making the proof of the Hessian 
case in Proposition A.l slightly asymmetric. 

Lemma A.4. Let f e C°(f2). If for every compactly-supported, smooth vector field W with \\W\\ e 
(7^(0), we have that f is a viscosity supersolution of the equation —W l W^ Hessy / + AijW l W^ = 0, 
then the function / satisfies Hess f < A in viscosity sense in the whole SI. 

Viceversa, if f e C° (SI) satisfies Hess / < A in viscosity sense in SI, then f is a viscosity supersolution 
of the equation —V % V^ Hessy / + AijV % V^ = Ofor every compactly-supported, smooth vector field V. 

Proof. Let us take a point x e Q, and a C 2 -function h in a neighborhood U of the point x such 
that h(x) = f(x) and h < f. Choosing a unit vector W x and a smooth, nonnegative function ip, 
which is 1 at X and zero outside a small ball inside U, we consider the smooth vector field W(y) = 
W x ip 2 (y), for every y € fi, which clearly satisfies || W|| = <p € C^(il). By the hypothesis of the 
first statement, the function / is then a viscosity supersolution of the equation — W l W^ Hessy / + 
A i jW i W j = which implies that -W*W£ Hess^ h{x) + Aij(x)W x Wl > 0. Since this holds for 
every point x € fl and unit vector W x , we conclude that Hess h(x) < A(x) as 2-forms, hence 
Hess / < A in viscosity sense in Q. 

The argument to show the second statement is analogous: given a compactly-supported, 
smooth vector field V, a point x E Q and a function h as above, the hypothesis implies that 
- V* VI Hcssy h(x) + Aij (x)V* Vi > 0, hence the thesis. □ 

Suppose now that / € C°(f2) satisfies Hess / < A in viscosity sense on the whole SI; hence, 
by this lemma, for every compactly-supported, smooth vector field V, the function / is a vis- 
cosity supersolution of the equation —V i yi HesSy / + AijV 1 V^ = 0. By Theorem A. 2 and the 
subsequent discussion, it is then a distributional supersolution of the same equation, that is, 



-fV%(V l V^)+A l3 V l V^ 



dx>0 



for every nonnegative, smooth function ip E C£°(f2). 

Considering a nonnegative, smooth function ip <G (SI) such that it is one on the support of the 
vector field V we conclude 

' /V;,il 'I r i./.r • / AijV^dx, 



ON THE DISTRIBUTIONAL HESSIAN OF THE DISTANCE FUNCTION 
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which means that Hess / < A in distributional sense. 

Conversely, if / G C°(fl) satisfies Hess / < A in distributional sense, then for every compactly- 
supported, smooth vector field W with \\W\\ C C^(fl) and smooth, nonnegative function ip € 
C£°(f2), we define the smooth, nonnegative functions ip n = (p + ip/n, where ip is a smooth, non- 
negative and compactly-supported function such that tp — 1 on the support of W. It follows 
that the vector field V = W^/ip^ is smooth, hence, applying the definition of Hess / < A in 
distributional sense, we get 



-fVUW l W>ip n ) + A V} W l W\ n 



dx>0. 



As (fin — > (p in (f2) and / is continuous, we can pass to the limit in n —> oo and conclude that 



/V;,iiril V; + .1,,U 'U'v 



> 0. 



for every nonnegative, smooth function ip e C£°(f2) and every compactly-supported, smooth 

vector field W with || W|| € C],(fl). That is, for any vector field W as above, we have that / is a 

distributional supersolution of the equation — W i W^ Hess^ / + A i jW l W : > = 0. 

By Theorem A.3 and the subsequent discussion, it is then a viscosity supersolution of the same 

equation and, by Lemma A.4, we conclude that the function / satisfies Hess / < A in viscosity 

sense. 

Summarizing, we have the following sharp relations among the weak notions of the partial 
differential inequalities Hess f < A and A/ < a, 

barrier sense ==> viscosity sense distributional sense. 
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